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Baryogenesis by heavy-neutrino decay and sphaleron reprocessing of both baryon 
and lepton number is reconsidered, paying special attention to the flavour structure 
of the general evolution equations and developing an approximate but sufficiently 
accurate analytic solution to the prototype evolution equation. Two different models 
of neutrino masses are examined, based on an Abelian U(l) or a non-Abelian U(2) 
family symmetry. We show that a consistent picture of baryogenesis can emerge in 
both cases, although with significant differences. 



Introduction to the revised version 

This paper contained two different original results: 

1 We discussed how to include flavour effects in computations of thermal leptogenesis, showing 
that one must write Boltzmann equations for a 3 x 3 density matrix that describes how the 
lepton asymmetry is shared between the 3 flavours. 

2 We gave approximated semi-analytical solutions to the standard Boltzmann equations and pre- 
sented numerical results in a way that covers all the parameter space systematically. 

Part 1 remains unchanged. Part 2 has been revised using the corrected Boltzmann equations of PQ. 
We do not list the changes we made. Final expressions are now simpler. 

1 Introduction 

The need for a mechanism to generate the baryon asymmetry of the universe after an early era of 
cosmological inflation is one of the main indications for physics beyond the Standard Model (SM). 
Several such mechanisms have been proposed, such as baryon-number violating interactions at tree 
level in Grand Unified Theories (GUTs), or sphaleron transitions taking place at the quantum level 
in electroweak baryogenesis. 

An even stronger indication for physics beyond the SM is provided nowadays by the anomalies in 
neutrino physics [211211], which can be interpreted in terms of neutrino oscillations. Using a simple- 
minded see-saw formula, m v = (200GeV) 2 /M, a neutrino mass between 0.03 eV and a few eVs, as 
seemingly implied by a coherent interpretation of the various neutrino results, suggests a mass scale 
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for new physics between 10 and 10 GeV. In turn, this mass scale could be related to the lepton 
number violating exchange of heavy right-handed Majorana neutrinos N. 

These recent findings in neutrino physics have stimulated the reconsideration of the possibility that 
the baryon asymmetry may be generated by a lepton asymmetry, arising from the out-of-equilibrium 
decay of the heavy N [£]. This conversion of asymmetries would occur by the reprocessing of both 
baryon (B) and lepton number (L) in sphaleron transitions. Ideally, from a detailed model of neutrino 
masses, both light and heavy, one would like to compute the difference between the number density 
of baryons and antibaryons, normalized to the entropy density of the universe, (rig — ng)/s = Yg, at 
the time of nucleosynthesis, known to be in the 10~( 10 ~ n ) range from the primordial abundances of 
the light elements In order to minimize uncertainties, the reheating temperature after inflation 
T r h is assumed to be bigger than the mass of the decaying neutrino [7j . 

With this general programme in mind, in this paper we add a few elements to the standard 
analysis of baryogenesis via leptogenesis: the consideration of the flavour structure of the problem 
and an approximate but sufficiently accurate analytic solution of the relevant evolution equations. 
These results would be fully relevant if a sufficiently detailed model of neutrino masses existed, which 
is not the case at present. Nevertheless, we apply our considerations to two different models of 
neutrino masses, based on an Abelian U(l) or a non-Abelian U(2) family symmetry respectively. A 
consistent picture of baryogenesis can emerge in both cases, although with significant differences, thus 
strengthening the view that makes leptogenesis an appealing mechanism for baryogenesis. 

2 General setting of the problem 

Sphaleron transitions, in equilibrium at temperatures above about 100 GeV, violate B and L while 
conserving the quantities Aj = \B — Li, i = e,/j,,r ITT] - As a consequence, by imposing the 
equilibrium conditions in the SM on all chemical potentials, one obtains below the electroweak phase 
transition [S] 

^ = pE y Ailo> Y Ai = ^Y B -Y Li , (2.1) 

i 

which specifies Yg in terms of the asymmetries YaJo> either initial or generated before the electroweak 
phase transition. 

In the pure SM Aj are conserved. Similarly J2i ^A* cannot be generated in a minimal SU(5) model. 
In order to generate a baryon asymmetry, we have to assume that an out-of-equilibrium interaction 
around a temperature T* violates (some of) the Aj and produces a net YaJo different from zero: N 
decay is the example that is of interest here, in which case T* is of order of the N mass, Mi. This 
decay acts as a source of asymmetry for the density of the lepton doublets Y^, related to Yj Ji by 
Yi i = Yi- + Y ei , where Y ei is the asymmetry of the right-handed charged leptons. If we only consider 
the ^-violating interactions, the Boltzmann equations for the have the (linearized) form 

Y ti = S i -7 i V ei , (2.2) 

where Yg. = dY^/dz, z = T/Mi, and both the sources Si and the wash-out coefficients 7$ are functions 
of z, properly normalized to account for the universe expansion. 

The way the evolution equations 1)2.2(1 are converted into evolution equations for the asymmetries 
Y\; El, defined in (|2.1|) . depends on which interactions are fast at T*, the decay temperature of N. 
Three temperature intervals are of interest: 

i) T J> 10 11 " 12 GeV, where the gauge interactions are in equilibrium while no lepton Yukawa cou- 
pling mediates equilibrium interactions. 

ii) 10 9 GeV <^ T <^ 10 11 " 12 GeV, where gauge interactions, the interactions due to the t, b, c, r- Yukawa 
couplings and (presumably) the sphaleron interactions are in equilibrium [lOj 1 . 

1 The rate of sphaleron interactions at finite temperature is uncertain. More recent works indicate an equilibrium 
temperature close to 10 10 GeV, lower than the value adopted here. If this is confirmed, it is straightforward to modify 
the following considerations. 
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iii) T < 10 9 GeV, where also the \x and s-Yukawa couplings and the CKM mixings mediate equilib- 
rium rates. 

The electron Yukawa interactions can be neglected as they are in equilibrium only at very low tempera- 
tures (T < 30 TeV) |H] • If we neglect the Y^.-violating interactions in eq. (|2.2|) , the different temperature 
intervals are characterized by different conservation laws, Yq u = 0. At the highest temperature (case 
i) the (globally) conserved quantities, other than Y^, include Yb, Yf*. and Y ei , which are progressively 
broken at lower temperatures. The different fast rates in any temperature interval give rise to dif- 
ferent equilibrium conditions for the asymmetries in the number densities, or the chemical potentials 

|H1 EU- For example fj>g T — [i er + fj>h = in the temperature range where the r Yukawa coupling 
\ T £ T e T h mediates a fast rate. These equilibrium conditions allow to express all the asymmetries in 
the number densities in terms of the Yq u , with expressions which depend on the temperature interval. 

It is now easy to see how the evolution equations 1)2.2(1 in presence of lepton number violation 
are converted into evolution equations for the Y&., directly relevant to the evaluation of the baryon 
asymmetry through (j2. 1|) . Since the fast interactions conserve Aj, they do not contribute to the 
equation for Y^., which takes the form 

^ = -($-7^), (2-3) 
where Y^ have to be expressed in terms of the Ya ; themselves, Yg t = Aij{T)Y^ j . Using the equilibrium 
conditions relevant to the different temperature intervals, the matrices of constant coefficients are given 
by 

i) A(T> 10 n - 12 GeV) = - diag(l, 1, 1); (2.4a) 



^(10 9 GeV^T^10 n ^ 12 GeV) = 34/351 -317/351 20/351 ; (2.4b) 



m 



diag(l,l,l) 
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A(T<10 9 GeV) = 29/506 -493/759 13/759 . (2.4c 



If the decay of the heavy neutrino occurs at the border of any of these temperature intervals, a more 
accurate treatment of the evolution equations would be needed. 



3 Flavour properties of leptogenesis by heavy neutrino decay 

The model of neutrino masses is specified by the Yukawa couplings of the heavy right-handed neutrinos 
N a , a = 1, 2, 3, to the left-handed SU(2)-doublets £i, i = e, fj,, r, and to the Higgs field h 

jO y = XiJiN a h + h.c. (3.1) 

which one can choose to write in the physical flavour basis of both the charged leptons and the heavy 
neutrinos of mass M a . In this flavour basis, the mass matrix of the light neutrinos is (v = (h) = 
174 GeV) 

v 2 

When needed, we accept here the currently dominant view that the light neutrino masses are hier- 
archical and that the two signals interpretable as due to the oscillations of three neutrinos are those 
related to the atmospheric and solar anomalies. In this context, the likely ranges for the eigenvalues 
of (|3.2j) , m„ 3 > m V2 > m Vl , are 

m vz = (AmL)' /2 ~ (0.03 - 0.1) eV; m U2 = (AmL) 1/2 < 10~ 2 eV (3.3) 

The flavour structure of (|3.1|) and (|3.2() is such that also the eq.s (|2.2() for Yg. are different in the various 
ranges of T* (~ Mi). As in the connection between Y^ and Y^, what counts are the temperature 
ranges described in the previous section: we discuss them in turn. For the time being, we consider 
the asymmetry produced by the decay of the lightest right-handed neutrino iVi . 
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3.1 T* « Mi > 10 GeV 

In this case the Yukawa and the sphaleron interactions can be approximately neglected during Ni- 
decay. Disregarding first the AL = 2 interactions, the evolution equation for Yi has the usual form 
discussed in the literature 

Y t = S-iYu (3-4) 
where S is proportional to the total decay asymmetry of N\ 

and f(x) ~ — 3/2x for x 3> 1 |12j . As long as the Yukawa interactions are negligible the lepton 
asymmetry can be considered concentrated on the £\ state, the lepton doublet to which N\ couples in 
the tree-level approximation: 

h = Xt&fyfiXXhi. (3.6) 

In an arbitrary flavour basis, the lepton asymmetry is described by a 3 x 3 matrix p arising from the 
difference of the density matrices of leptons and anti-leptons and normalized so that Tr p = J2i Yi t (13] . 
In the case we are considering, as shown below, p is proportional to the projector Pi on the state £\ 
up to non-diagonal contributions. When the Yukawa interactions come into equilibrium at lower 
temperatures, they simply kill the off-diagonal terms of the matrix p in the physical lepton- flavour 
basis, while leaving unaltered the trace, which is what influences Yli or > ultimately, the baryon 
asymmetry. 

Let us now come to discuss the possible washing effect of the AL = 2 interactions mediated by 
off-shell iV a -exchanges. Their rates, at T < Mi, are controlled by the light neutrino mass matrix (|3.2|) . 
Focusing on the heaviest eigenvalue, ()3.3|) . the dominant AL = 2 amplitude will act on the state v-$ 
and will be proportional to m Us /v 2 . The corresponding rate is in equilibrium at temperatures above 
10 n ~ 12 GeV, so that it is only in this range that the lepton asymmetry produced by Ni is affected. 
This effect is readily included in the evolution equation for Yg, eq. (|3.4|) . if £3 coincides with £\, eq. ()3.6|) . 
as it happens if the exchange of N± dominates the mass matrix (|3.2|) of the light neutrinos. This is 
usually done in the literature Jl] by introducing a properly normalized AL = 2 rate jal=2 in the 
left-hand-side of eq. (|3.4jl . which becomes 

Y l = S-( 1 + lAL=2 )Y l . (3.7) 

The inclusion of the AL = 2 washing effect is different if £3 is not aligned with £±. If we keep only 
the dominant AL = 2 washing interactions acting on £3, their effect is accounted for, in the general 
case, by projecting the evolution equation for the matrix p in the subspace £3 and in its orthogonal 
complement, and writing ^Y^ = Y3 + Yj_ with 

Y 3 = c 3 S- [Tr(PiP 3 )7 + 7AL= 2 ]>3 (3.8a) 
Y± = (l-c 3 )5-7[l-Tr(PiP 3 )]Y ± . (3.8b) 

Here Y± describes the trace of the matrix p restricted to the corresponding subspace, 

_ r(iVi -+ £ 3 h) - r(iVi -> i 3 h) 
03 r(JVi -> eh) - r(iVi -+ih) [ ' ) 

is the part of the asymmetry produced by the decay of Ni into £3 and Pi, P3 are the projectors over 
£1, £3 respectively. 

This intuitive result can be formally derived by writing the evolution equation for the matrix p 
associated with the lepton asymmetry 

p = S{P l + — )-7 — 7al=2 — g — , (3.10) 
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where AP/2e\ in the source term accounts for the misalignment in flavour space of the states £ and 
£ to which iVi decays, with projectors P and P (see appendixEJ). An exact expression for the source 
term is 

T(N^fh)P-T(N^£h)P 
T(N -» £h) - T(N -» £h) ' y ' ' 

which, expanded to first order in the asymmetry, gives the source term in (|3. lOf) with AP = P — P 
and ei given by eq. (|3.5j) . It is important to notice that AP/2e\ is of order unity and cannot be 
neglected. Without 7al=2, since \{P\, AP} = AP + 0(e\ ), the solution of (|3.1U|) to first order in the 
asymmetry is 

"= r '( Fl+ fr) (3 - 12) 

where Y t satisfies (JOJ. Notice that Trp = Y e and that Tr[Pi(Pi+ AP/2ei)] = l + 0{e{). With 7al=2, 
the restriction of eq. (|3.1U|) to the two orthogonal spaces done before leads to eq.s ()3.8aj) and 1)3. 8b|) 
since (see appendix |A")) 

Tr\pJ Pl + ^l 



c 3 . (3.13) 



3.2 10 9 GeV < T* & M x < 10 11 - 12 GeV 

In this case, during Aq-decay, the AL = 2 interactions can be neglected. On the contrary, one has to 
take into account the fast rates due to the r- Yukawa interactions. One of their effects is best seen in 
the flavour basis, where they drive to zero the off-diagonal terms p Te , p Til of the matrix p. Defining 
c T in analogy with C3, eq. ()3.9|) . one would have 

Ptt = Yi T = c T S - 7 Tr(P 1 P T )Y tT (3.14) 

and, for the sum of the electron and muon asymmetries, indistinguishable at this stage, 

Pw + Pee = % + Y ie = (1 - c T )S - 7 [1 - MPlPrWl, + *U (3.15) 

Furthermore, the equilibrium of the r- Yukawa and sphaleron interactions has the effect described in 
section|21and accounted for by the 2x2 matrix A nm which mixes A T with A e + A^ (re, m = {e + p, r}). 
In analogy with section |21 

/(34 - 317)/351 40/351 
Anm ~{ 1/117 -82/117 



(3.16) 



3.3 T*^M!<10 9 GeV 

The extension to this case, when both the r- Yukawa and //-Yukawa interactions are in equilibrium, is 
obvious. All the off-diagonal elements of the matrix p in the flavour basis are driven to zero so that 

Y e .=c i S- 1 Tr(P 1 P i )Y ii (3.17) 

where i = e,p,T. Including the fast Yukawa and sphaleron interactions as described in section |21 with 
the proper "mixing" matrix ()2.4b|) one obtains the evolution equations for Y^. 



4 Approximate analytic solutions of the Boltzmann equations 

In this section, and in the related appendix El we describe an approximate but sufficiently accurate 
analytic solution to Boltzmann equations of the form (|3.7j) . Numerical solutions have been presented 
in where the relevant SM rates have been computed. As previously discussed, this equation is 
appropriate for T* > 10 11 " 12 GeV and when the exchange of iVi dominates the mass matrix of the light 
neutrinos. Boltzmann equations with a similar structure hold at lower temperatures, when sphaleron 
and Yukawa interactions are fast. Similar analytic solutions can be developed in these cases. 
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Figure 1: Contour plot of the efficiency factor r/ of leptogenesis in the SM as function of 
(fhi,Mi). The AL = 2 interactions not mediated by N\ have been computed assuming m U3 = 
max(mi, (3r 10- 3 ) 1 / 2 eV) and assuming that their flavour structure gives the weakest washing (X = 1) 
in fig. Ufa, and the strongest washing (X = m 2 /mfj in fig-US). Continuous (dashed) lines assume ther- 
mal (zero) initial abundance of N\. 



In the most general case, N\ exchange gives a contribution fh\\(l\h* /v) 2 +h.c] to the light neutrino 
mass operator 2 : the flavour of l\ can be different from the flavour of the heaviest neutrino vj, in the 
is doublet, and fh\ = (\^\)nv 2 /M\ can be smaller than its mass m„ 3 . The appropriate system of 
eq.s Q3.8JI depends on the relative flavour between l\ and £3. They simplify to a single equation of the 
form (|3,7|) in the two extreme limits: 

(a) l\ and £3 have 'orthogonal' flavours (in the notation of section |3] this means Ti(PiP^) = and 
c 3 = 0); 

(b) l\ and £3 have 'aligned' flavours (Pi = P3 and C3 = 1). 

In case (a) the dominant AL = 2 interactions cannot wash out the leptonic asymmetry in £\, as would 
happen for example if they only acted on \i and r flavours while t\ = £ e . The relevant equation ()3.8b|) 
has the form ()3.7|) . with jal=2 given by the subdominant AL = 2 interactions mediated by N\ 
(neglected in (|3.8|0 . In case (b) the relevant equation (|3.8a|) has the form (|3.7j) . with the AL = 2 rate 
7AL=2 being X = (m„ 3 /mi) 2 times stronger than the one mediated by N\ alone. The intermediate 
situation can be studied by solving (|3.8|) for any given relative flavour misalignement between £\ and 
£3: the result should be intermediate between the ones obtained for X = 1 (case a) and X = m^/fhf 
(case b). 

Both the source 5 and the decay coefficient 7 + 7al=2 in eq. (|3.7|) depend upon , the N\ 
density relative to the total entropy density, which satisfies its own evolution equation. We define (see 
appendix lB|) 

ml = 2.3 10" 3 eV. (4.1) 

2 Our parameter mi is renormalized at the high energy scale ~ Mi. High energy neutrino masses are r ~ 1.2 1.3 
times larger than at low energy 1241 IT) . 
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If fh\ -C m\, the decay width of N\ is much smaller than the Hubble constant and N\ decays strongly 
out of equilibrium, while, conversely, it is almost in equilibrium for fh\ S> fh\. Explicit expressions 
for Yn-i(z) in the two regimes are given in appendix IE! 

With explicit knowledge of S and 7, eq. (|3,7j) can be integrated to obtain Yi at T <C M\ 



Y e (z ~ 00) = e iY£(0)u i.e. M = ^ oo) = 



where 77 < 1 is an efficiency factor describing the effect of wash-out interactions. In appendix [B] 
we describe an accurate analytic approximation for 77. The first-order linear differential equation for 
Yi, eq. 1)3. 7 Jl . can be explicitly solved in terms of integrals. Our approximation for 77 is obtained by 
inserting the asymptotic expressions for the interaction rates and for the N\ abundance described in 
appendix IB . 1 1 and evaluating the integral in the saddle-point approximation. 

Depending on the mass and Yukawa couplings of N±, four different regimes are possible. The 
AL = 2 interactions can be either significant or negligible, and the N± decay rate (which is comparable 
to the thermally averaged AL = 1 interactions) can be faster or slower than the expansion of the 
universe, depending on the value of fh\/fh\. Specializing our general approximation ()B.5j) described 
in the appendix to the four regimes, we can derive less precise but more explicit expressions for 77: 

A. N\ decays strongly out of equilibrium and all the wash-out interactions are negligible. In this 
case 77 1. 

B. N\ decays almost in equilibrium and AL = 1 interactions are not negligible. As explained in 
appendix El 77 is suppressed only almost linearly in fh\/fh\\ 

77 w 0.168 11 — - — (4.3) 



We put a numerical coefficient 10% lower than what suggested by HB.5|) so that this approx- 
imation agrees with our full computation within few %. Notice that 77 depends, up to small 
corrections, on the mass and Yukawa couplings of N\ only through the single combination fh\. 

C. N\ decays almost in equilibrium and AL = 2 interactions are not negligible. These interactions 
wash out Yi exponentially: 



77 ~ exp 



mi / M\X 



m* V 3.3 10 15 GeV 



(4.4) 



The factor in the exponent is proportional to the square root of 7al=2 and gives a significant 
suppression only if Mi is heavy enough. 

D. Ni decays strongly out of equilibrium and AL = 2 washing is not negligible. These interactions 
again wash out Yi exponentially: 



77 « exp 



m) V0.8 10 13 GeV ml /mj 



(4.5) 



If X = 1, this suppression is significant only for uninteresting high values of M\. 



These four regions are clearly visible, and explicitly indicated, in fig. where we show contour plots 
of 77, as function of the two unknown parameters 771-1 and Mi in the two extreme situations alluded to 
above: fig. [^i refers to case (a) and fig. to case (b). 

Before concluding this section, we make a few comments on the validity of our approximations. 

1. We have assumed that Ni has a thermal abundance at T S> Mi but this could not be the case. If 
7Tii S> m\, Ni rapidly thermalizes anyhow, so that its 'initial' abundance is irrelevant. If instead 
7Tii <C m*, thermalization is slow and leptogenesis becomes sensitive to the initial abundance: 
starting with Y/y 1 (0) = 0, the Ni abundance at decay is roughly suppressed by a factor rhi/m\, 
not included in fig.s ^ 
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2. In presence of a strong exponential suppression, the slow (i.e. 7 < 1) Yukawa and/or sphaleron 
interactions can not be neglected: a fraction rj ~ 7/7^ of the asymmetry in left-handed leptons 
I is transmitted to right-handed leptons and/or to quarks before being washed out. 

3. We have done the computation in the SM. Apart from factors of order one, there should be 
no significant difference between the SM Jl] and MSSM ^2] predictions for Yg, with the main 
change possibly due to a large value of tan (3. 

4. We have neglected the asymmetry generated by the heavier neutrinos, even if their mass is below 
T r h- This asymmetry is generally believed to be washed out by fast interactions mediated by 
the lightest right-handed neutrino. Considering the flavour structure of the problem we can see 
that this is not always the case. The heavier neutrinos produce an asymmetry of the generic 
form (see eq. (|3.12|) ) 

» = li(ft + ^). (4.6) 

For T* > 10 11 " 12 GeV, analogously to eq. (|3.8aj) and (|3.8b|) . only the restriction of ph to the 
subspace l\ will be washed out by N\ interactions, while its orthogonal complement will remain 
untouched. For lower T* the wash-out efficiency depends on the projection of and i\ on the 
flavour states: also in this case a non-negligible part of the asymmetry can survive. 



5 Baryogenesis in specific models of neutrino masses 

As made clear from the previous discussion, a detailed calculation of the baryon asymmetry generated 
via leptogenesis would require a detailed knowledge of the model for neutrino masses, not available at 
present. A typical model is currently specified by textures for the matrix A of the Yukawa couplings and 
for the mass matrix M of the heavy right-handed neutrinos, with undetermined complex coefficients 
of order unity for every non-zero entry. The ignorance of these coefficients induces an uncertainty in 
the lepton asymmetry, generically of the same order unity. A pile up of effects correcting the final 
result for Yi by up to one order of magnitude in either direction, although not likely, cannot be safely 
excluded. 

Before considering specific models, it is useful to recast the relevant equations in a more expressive 
form, emphasizing the role of the relevant parameters [TJ1|- From ()3.5j) . the total decay asymmetry 
can be written as 

N = T7, -?—9 5.1 

l07T V z 

where g is a model dependent flavour factor which, barring cancellations, is less than unity and does 
not depend on any absolute scale. Ref. J7j shows that g < 1 if Mi <C M2,3- We find therefore, 
from (fl?2|) 

IO-10 ~W i 9 GeV [ ' 

where the efficiency factor 77 is also less than unity and depends on the coefficient fh\ of the contribution 
mi[(£ih* /v) 2 + h.c] to the mass matrix of the light neutrinos from N\ exchange. 

In the light of the above general considerations, we compute the expected baryon asymmetry in 
two models of neutrino masses, based on an abelian U(l) or a non-abelian U(2) family symmetry 
respectively. In both cases we have checked that the careful treatment of flavour, as described in 
section |21 and |2J can modify the result with respect to the naive treatment based on eq. (j2.2j) only to 
an amount which stays within the model uncertainty. Hence we only describe the results of the naive 
analysis. As in all the previous discussion we assume a reheating temperature after inflation above 
Mi. 
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5.1 The U(l) model 



The essential elements of this model JB] are the U(l)-charges of the three left-handed lepton doublets 
^3> £\'- a, a, a + x, and the biggest charge 0\ of the heavy neutrinos, associated with iVf. The tilde 
on top of I and N c is there to remind that they are not mass eigenstates. The U(l) charges of all 
particles with the same chirality — hence the 1C ' on the heavy neutrino fields — are non-negative and 
are taken consistent with SU(5) unification, so that they are identical within each SU(5) multiplet. 
Every element Aj a of the Yukawa matrix is given, up to a factor of order unity, by S q ' +ea , where 
5 is a small parameter and qi, a are the U(l) charges of the fields £i and N° respectively. An 
analogous formula holds for all other Yukawa couplings. In the same way, again up to factors of order 
unity, the mass matrix of the heavy neutrinos is M a p ~ Mq S^ a+9l3 \ where Mo is an unspecified 
mass scale. The equality of the U(l)-charges of £3 and £2 is at the origin of the large neutrino 
mixing supposedly observed at SuperKamiokande. The large top Yukawa coupling and consistency 
with SU(5) unification require the vanishing of the U(l) charge of the third generation right-handed 
charged leptons, approximately tr. 

With these inputs, from the explicit construction of the A and M matrices, it is immediate to 
obtain 

si^^S 2 ^, m„«mi«^i M^S^Mo. (5.3) 

1D7T Mq 

Phases of order unity are assumed in the various matrix elements. We have also assumed that the 
unspecified charges of the heavy neutrinos other than Ni are sufficiently smaller than 9\ so that 
M~i/M a ~ (5 2 ( e i- a) 5 a. = 2,3, are small. A posteriori, this can only be marginally the case. A 
drawback of the model is the need to invoke a fortuitous cancellation to explain the hierarchy in 
the neutrino mass matrix, m U2 /m U3 <0.05, against the naive expectation m U2 « m„ 3 . In its unified 
version, the model also gives m T ~ 5 a v, although with an arbitrary identification of the vacuum 
expectation values of the Higgs fields relevant to m T and m V3 . 

Eq.s (|5.3|) can be combined to give the flavour factor g in (|5.1|) . approximately equal to unity. 
Furthermore Mi, using also m T ~ 5 a v, is bounded by M 1 < M « m 2 T /m U3 < 10 11 GeV. There 
is no other restriction on M\ which can therefore be fixed to obtain the desired value of Yb from 
fig. ([T]). Since rh\ ~ m U3 ~ (0.3 -j- 1)10 _1 eV, from fig.s^ rj = (0.3 -r- 1)10~ 2 . Hence consistency with 
Y B = (10~ 10 ^- 10" 11 ) is obtained for Mi w lO 10 ^" 11 GeV. Although lower than M < 10 11 GeV, this 
value is such that also the heaviest singlet neutrinos have to be close in mass and all in the 10 11 GeV 
range. 



5.2 The U(2) model 

The relevant pieces of information are contained in the form of the matrices for the £N Yukawa 
couplings and for the heavy neutrino masses 



1 



e e e 
e 1 



M « eM 







e' 1 1 
1 e 



(5.4) 



where e ~ 0.02 and e' ~ 0.004 are dimensionless parameters related to the hierarchical breaking of 
U(2) and determined from the observed values of charged fermion masses |19j . The tilde on top of A 
indicates that A, unlike A in eq. (|3.2|) . is not in the physical basis for the heavy neutrinos but rather 
in the basis specified by M. As before, order one prefactors are understood in every entry of (|5.4|) . 
whereas the blanks stand for non-zero but sufficiently small elements. Since the structure of the model 
is intimately tied to unification, eM is of order M GU T ~ 2 • 10 16 GeV fTHlEHj . 

The diagonalization of M in eq. (|5,4|) leads to M2 ~ M3 ~ cMq ~ Mqut and to a substantially 
lighter Mi ~ ee /2 MGUT ~ 10 9 -v- 10 10 GeV. In turn the £-N Yukawa couplings in the physical N-basis 
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(N±, N2, N3) acquire the form 



£ /2 £ / £ / 



A « I e' e e I . (5.5) 
e' 1 1 



The exchange of M\ dominates the light neutrino mass matrix, giving 



0.1 eV, 1^ « e « 2 ■ HT 2 , (5.6) 



M G ut e m 



V3 



in reasonable agreement with observations. Furthermore mi « , so that 77 « (0.34- 1)10~ 2 and, for 
the flavour factor in (|5.2|) . 5 « m V2 /m U3 « e « 0.02. Putting everything together, we get from (|5.2I) 



10 8 GeV 



Unlike the previous case, however, Mi is not a free parameter. Since Mi » 10 9 - 10 GeV, the baryon 
asymmetry is fixed to 3 10~( 12H ~ 13 ). 



6 Conclusions 

The experimental results in neutrino physics of the last few years have made plausible the case for a 
baryon asymmetry generated by a lepton asymmetry arising from the out-of-equilibrium decay of a 
heavy right-handed neutrino. With this motivation, this paper achieves two purposes. The first is to 
determine in a precise way the evolution equations for the baryon and lepton asymmetries by taking 
into account the full flavour structure of the problem. This would be fully relevant with a sufficiently 
detailed model of neutrino masses on hand. 

The second purpose is to examine the expected asymmetry in two specific models of neutrino 
masses, within their limit of uncertainty In both cases, we have shown that a consistent picture of 
baryogenesis can emerge, although with significant differences. In the U(l) model, which does not 
account automatically for the hierarchy of neutrino masses, the observed baryon asymmetry can be 
obtained by choosing the scale of the heavy neutrinos, all close in mass, between 10 10 and 10 12 GeV. In 
the U(2) model, the lighter of the heavy neutrinos has a mass Mi fixed by independent considerations 
at Mi ~ 10 9 ~ 10 GeV and significantly lower than that of the two others M2 ~ M3 ~ Mq. In this case 
the baryon asymmetry is fixed at 3 10~( 12 ~ 13 ) against an observed value of lO - ^ 10 " 11 ). In view of the 
uncertainties discussed in the text, we consider this as a success, strengthening the view that makes 
leptogenesis an appealing mechanism for baryogenesis. In every instance it is clear that the decaying 
neutrino is heavy, above 10 9 GeV and correspondingly T T ^ is bigger than the same value, a non trivial 
constraint on the evolution of the early universe. 
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CT96-0090 and contract No. ERBFMBICT983132. A.S. thanks G. Giudice, A. Notari, M. Plumacher, 
M. Raidal and A. Riotto for useful discussions. 



A The explicit form of Ni decay amplitude 

The tree-level amplitudes and the one loop corrections to N\ decay are, up to an overall factor 

en = a(Ni -> tji) = Xn 5j = Aa(Aq -> tji) = \i a XjiXj a A a 

at ee a(Ni hh) = X* n 8 { = Aa(JVi iih) = \l a \ n \* ja A a , 
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where lmA a = -f(M a /M 1 )/32ir and / is defined in eq. Note that Si / 5*. Neglecting AL = 2 

interactions, the Boltzmann equations for the density matrices of leptons and antileptons are 

p e = S e P + j(P-^{P, Pe }) (A.la) 
n = S l P + 'y{P-±{P,pi}), (A.lb) 



where 

_ (a + 6)i(a + 5)] _ (a + 6)i(a + S)) 

13 ~ \a + 5\ 2 13 ~ \a + 6\ 2 ' { > 



By taking the difference between (jA.laJ) and (jA.lbj) . defining S = Si — Si and linearizing in p = pg — pg, 
eq. (|3.10|) follows. In eq. (|3.10|) we have also included the AL = 2 terms, which can be discussed along 
similar lines. Explicitly we have 

p APN i Im(AqA* a A^A fcl )A a 

1 2ejy 2EJm(A l iA* Q A^A fc i)W Q l ' j 

so that eq (|3,13|) is easily justified. 

B Approximate solutions of Boltzmann equations for leptogenesis 

As usual |211 1141 ITo] we assume a Boltzmann kinetic distribution for bosons and fermions: the Boltz- 
mann equations for the total abundances n^, ng and ng are expected to be correct within few 10% 
errors. With this approximation the Hubble constant at temperature T, as predicted by standard 
cosmology, is given by H 2 {T) = ?fG N p where M P1 = G N 1/2 = 1.22 10 19 GeV, p = 35 SM T 4 /vr 2 and 
gsM is the number of ultra-relativistic spin degrees of freedom (gsM = 118 in the SM at T 100 GeV). 
The entropy density is s = 4gsM^ 3 /tt 2 . The number density of a particle p with g p spin degrees of 
freedom (for example gN 1 = 2 and gg i = 2) and mass M p in thermal equilibrium at temperature T is 

„ w = e, (a) = ™ (0) _ _Jt_ 



where z = M p /T and K^z) is a Bessel function (precisely defined as BesselK[2,z] in Mathematica 
notation [221) with the following limiting behaviors: K 2 (z) — 2/z as z — ► and K 2 (z) ~ e~ z yf^j2z 
as z — ► oo. The final 1? asymmetry can be written as 

s 37 JVl 37c?sm 

where ei is the CP asymmetry in N± decays given in eq. 1)3.5(1 . and 77 is an efficiency factor for 
leptogenesis that depends on the mass and Yukawa couplings of N± and should be determined solving 
the relevant set of Boltzmann equations. We assume thermal equilibrium at T ^ M\: consequently 
rj < 1. In this section we derive the approximate analytical expression for rj used in section 0] As 
usual it is convenient to write equations for Y^ 1 (z) = n^^z) / 's(z) and for Yg = (ng — where 
z = M\jT and s is the total entropy density. 

B.l The Boltzmann equation for the Aq abundance 

We neglect corrections of relative order C(A 2 , g 2 )/^ 2 (some of them are included in state of the art 
computations). A first equation 

Y^+ lNl (Y Nl -Y^) = where OP 7*1 = ( B1 ) 

szn 1 Ni 
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controls how much out of equilibrium Ni decays. All quantities depend on z = M\/T. An approximate 
expression for the dimensionful thermally averaged decay rate jd can be found at the end of this 
appendix. 

Assuming thermal equilibrium of N\ above the decoupling temperature T ~ Mi, the N\ abundance 
at decoupling can be either almost in equilibrium (if 77^ S> 1) or strongly out of equilibrium (if 
77V t <C 1). Approximately 7^ ~ mi/mf, where 



1 M Pl (l+g Nl ) 

depends on cosmology (the numerical factors are appropriate for later use). All the dependence on 
the mass and Yukawa couplings of Ni is incorporated in mi = Aft> 2 /Mi, the contribution to the 
light neutrino mass mediated by N\ (A 2 = (A'A)m is the 'total' squared Yukawa coupling of Ni). 
Approximate solutions for valid in the two regimes are: 

• If fh\ <C mf , Ni decays strongly out of equilibrium. In this case the only relevant collision term 
is the thermally averaged decay rate at T <C Mi, so that the Boltzmann equation for N\ and its 
solution are 



. <y D (z ~ 00) Y Nl v c q^n^ 
Nl = 7z~H y nAz) =Y N ^(0) exp 



mi 4v2/vr 



mi 



At very small rh\ N\ decays slowly at T^m^/gsM an d gives a substantial contribution to the 
total energy density: we do not study this possibility PQ. 

If fh\ 3> fh*, Ni decays almost in equilibrium. If the collision terms 7^ are fast, the approximate 
equation for A(z) = [Y^iz) /Y^{z)} — 1 and its solution are [U 

Y /eq H 



with Y*? = -z 2 Ki(z)/4gsu- 



Even if we do not have a unique good approximation for the N\ abundance, we can obtain a good 
approximation for the generated lepton asymmetry. The reason is that, if the washing interactions are 
effective (otherwise 77 is trivially close to 1), only the lepton asymmetry produced by late N\ decays 
survives. The final result is insensitive to what happens in the range of temperatures T ~ Mi, where 
we do not have a simple approximation for Y^ 1 (z). 



B.2 The Boltzmann equation for the leptonic asymmetry 

We now consider the Boltzmann equations for the leptonic asymmetry Yi = (ri£ — n^)/s. For concrete- 
ness, we give approximate solutions to the equations usually employed in the literature. If sphalerons 
and/or lepton Yukawa couplings and/or a non minimal flavour structure give significant interactions 
at T « Mi, the equations must be modified by appropriate order one coefficients. The Boltzmann 
equation for Yi is 

Yl + ll {z)Y l = Si{z), (B.2) 
where the source and damping factors Si and ji are given by 

St = ^, 7 . = 7D + 47 g- (B.3) 
1 szH ' U 2szHY c cq V ' 



12 



Approximate expressions for these factors, inserting our almost-in-equilibrium approximation for 
Yn 1 (z) (it is immediate to pass to the other case), are: 



S t (z) « 1 - £l Y^^e l ^z^e- z (B.4a) 



4gsM 

z>l 2m 

3^* Z " 6 +_ z 2 7T 5 /: 



^ 5/2e - z + 3V2Ai x 



(B.4b) 



The term in ()B.4b|) decoupling as e~ z at low T is due to AL = 1 interactions, while the term multiplied 
by X is due to AL = 2 interactions and decouples only as 1/z 2 . As explained in section 31 the factor 
X, ranging between 1 and m 2 3 /m 2 , takes into account possible significant AL = 2 interactions not 
mediated by Aq. 

The Boltzmann equation for Yi is a linear first order differential equation: using its standard 
solution in terms of integrals we finally obtain 



roc 



r]=J o dz> expl-F(z')] where F{z) = - In + / 7^')^'- (B.5) 



Since e~ F is a bell-shaped function, by integrating the 2nd order Taylor expansion of F(z) around 
its minimum (F'(z) = 0) gives a more explicit but less accurate expression, r\ ~ e~ F ^ \J1-k / F" (z) . 
When AL = 2 scatterings are irrelevant this expression, with Si and ji taken from (|B.4|) . gives a good 
approximation even if N± decays strongly out of equilibrium. 

Depending on the mass and Yukawa couplings of N\, four different regimes are possible. Special- 
izing our approximation (|B.5|) to the four regimes, we have derived the less precise but more explicit 
approximations for r/ presented in the text. It could be a bit surprising that, when the 7al=i inter- 
actions are much faster than the expansion rate of the universe, the lepton asymmetry is suppressed 
only linearly. The reason is that AL = 1 interactions, having Aq as an external state, at low tem- 
perature are suppressed by a Boltzmann factor: ^/al=i(z) ~ 7AL=i(0)e~ z . For z>z these washing 
interactions become negligible and all the Aq decays give rise to unwashed leptonic asymmetry, so 
that the suppression factor is approximately given by 

77 « Y Nl (z)/Y Nl (0) « e~~ z « l/ 7 AL=i(0) ~ fhl/mv 

The AL = 2 interactions mediated by N{ can be computed in the effective theory below Aq and 
decouple only as 1/z 2 . These interactions do not decouple exponentially and consequently can wash 
out Yi exponentially. If fh\ ~ m Us these interactions are significant when Ai ~ 1. 

B.3 Approximate expressions for the thermally averaged interactions rates 

Evaluating the various thermally averaged interaction rates Jl] is the more difficult step of a numerical 
computation. Quite simple and accurate expressions hold in the low temperature limit 

MfXl e- sub (MjAi) 4 

where X 2 = (XX^u is the 'total' squared Yukawa coupling of A^, and we have followed the notation of 
for the dimensionful 7: -)d is the thermally averaged decay rate, 7jy b represents AL = 2 interactions 
mediated by N. As shown in [Q there is no quasi-resonant enhancement of 7|y b at T ~ mjVi so 
that we can neglect it. These approximations have been used to obtain eq.s (|B.4j) . In our analytic 
approximations we neglect AL = 1 scatterings, that at small temperature give corrections of relative 
order 0(X 2 t , ffD/vr 2 , e.g. 

3M 4 x2 , e~ z , , 3Mf l2x2 e" z , M x /Ti „, , 

jHs{z ~ 00) « — -^=A 1 A t - — 7Ht(z ~ 00) sa A x A f - — -y^ln B.6b 

16V2 {-Kzfl 1 64v 2vr 2 (7rz) 5 / 2 m/j 

AL = 1 scatterings are subleading also at large temperature T 3> M\, because inclusion of thermal 
effects modify 70 oc T 2 into 70 oc T 4 [I] . 
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